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ABSTRACT 


A  sea  floor  benchmark  experiment  was  conducted  in  an 
area  about  16  nmi  west  of  Pt.  Lobos,  California  (36°30'N  x 
122°17'  W)  during  18-22  May  1985. 

Two  baseline-crossing  methods  were  used  to  determine  the 
relative  positions  of  acoustic  bottom  transponders.  The 
method  of  least-squares  adjustment  was  used  to  analyze  the 
data.  Relative  position  determination  of  the  transponder 
array  is  discussed  and  recommendations  are  made  for  futher 
improvement.  The  advantage  of  these  methods  is  their 
simplicity.  Their  disadvantage  is  the  relatively  large 
amount  of  ship  time  they  require  to  achieve  acceptable  accu¬ 
racies. 

Transponder  arrays  such  as  the  one  deployed  can  be  used 
for  solving  many  types  of  problems  in  sea  floor  engineering, 
which  will  be  of  increasing  importance  in  the'future. 
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I.  IHTBODOCTIOH 


A.  BACKGfiOOHD 

In  today’s  world,  the  ocean  scientist  and  engineer  need 
increased  accuracy  for  marine  positioning  to  accomplish 
their  work-  The  establishment  of  a  system  for  the  precise 
determination  of  positions  on  the  deep  seafloor  is  one  of 
the  most  fundamental  challenges  of  working  in  the  open  sea. 

A  seafloor  benchmark  positioning  system  is  proposed 
which  can  aetermine  accurately  the  location  of-  objects  on 
the  sea  floor  (Saxena,  1974).  Such  precisely  located  beacn- 
marks  can  later  be  used  in  turn  to  delineate  offshore  prop¬ 
erty  lines  and  national  boundaries-  They  are  also  useful  in 
solving  problems  in  seafloor  engineering,  in  plate  tectonic 
studies  and  in  connection  with  bore-hole  reentry  associated 
with  off-shore  oil  recovery. 

The  objective  of  this  study  is  to  demonstrate  the  feasi¬ 
bility  of  establishing  such  seafloor  benchmarks. 

B.  BESCRIPTIOH 

Early  experiments  associated  witn  seafloor  acoustic 
transponder  arrays  were  carried  out  by  Hart  (1967),  Haehnie 
(19o7),  Fubara  and  Mourad  (  1972),  ana  JlcKeown  (1975),  wno 
used  a  naseline  crossing  method  to  solve  for  the  distances 
between  transponders.  "Baseline"  in  this  taesis  means  tne 
line  between  the  projections  of  two  bottom-mounted  acoustic 
transponders  onto  the  sea  surface- 

Ihe  baseline  accuracies  obtained  by  Hart  (1967),  acKeown 
(1975),  and  Fubara  and  Hourad  (1972)  were  ±15-7  m  for  base¬ 
lines  of  5509  m,  ±3.  1  to  ±4.2  m  over  baselines  of  6373  to 
7219  m  and  ±15.5  m  over  baselines  of  9364  m,  respectively. 


The  first  two  of  the  above  authors  did  not  use  a  least 
squares  method  for  data  analysis. 

There  are  three  parts  to  the  determination  of  the  abso¬ 
lute  position  of  an  ocean-bottom  transponder  array  using 
satellites.  They  are  (Figure  1-1): 

1.  To  determine  the  geodetic  position  of  tne  ship’s 
receiving  antenna  from  the  Global  Positioning  System 
(GPS)  ; 

2.  To  convert  the  geodetic  position  of  the  antenna  to 
tne  ship- mounted  acoustic  transducer;  and 

J-  To  determine  the  geodetic  position  of  the  acoustic 
transponder  array  on  the  sea  floor  using  a  snip- 
mounted  acoustic  transducer. 

This  thesis  is  concerned  with  the  relative  position  of  a 
bot tom- mounted  transponder  array  with  regard  to  the  ship- 

C.  DESCfilPTIOil  OF  THE  SYSTEM 

Acoustic  transponders  manufactured  by  Oceano  Instruments 
were  used  for  this  project.  They  were  deployed  by  the  S/V 
Acania.  Figure  1.2  shows  the  configuration  of  the  Oceano 
system  used-  Their  principal  features  are: 

1-  Ail-121  :  (Acoustic  Module) 

The  hydrophone  was  fixed  amidships  to  tee 
port  side  at  a  depth  2.95  a. 

-  Transmission  frequencies  :  8  to  16  k.Kz- 

-  Reception  frequencies  :  6  to  16  icHz. 

2.  RT-121  :  (Recoverable  Transponder) 

-  Frequencies  :  8  to  16  hHz- 

-  Delay  time  :  15-0t0.1  ms. 

TT-2Q1  :  (The  Accoustic  Telecommand  Module) 
has  the  following  functions; 

-  Coce  acoustic  signal  transmitted  by  either  the 
acoustic  module  or  UQC  transducer. 

-  Monitors  range  to  desired  transponder- 
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Figure  1.  1 

The  Configuration  of  Bottom  Transponder 
Shown  with  Ship  and  GPS  Satellite. 
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-  Monitors  receptioa  of  code  by  traaspoader. 

.  IM-100  :  (Data  Interface  Module) 

-  Accepts  up  to  16  data  inputs- 

-  Max-  converted  data  format  ;  7  digits  without 
sign  or  6  digits  plus  sign. 

-  Computer  interrupt  request  possible  through 
any  input. 

-  Each  channel  may  be  read  separately. 

5.  RM-201  :  (Rangemeter) 

-  Transponder  turnaround  delay  offset  :  0  to  99.99  ms. 

-  Mean  velocity  of  sound  input  :  0  to  9999-9  units  per 
second. 

-  Standard  receptioa  filters  :  8  to  16  k,Hz  (0.5  kHz  step) 

-  Bandwidth  :  203  Hz. 

Critical  to  the  experiment  are  the  following: 

1,  The  transducer  (a  part  of  the  AM-121)  was  2.95  m 
below  the  surface,  and  the  transmission  and  reception 
frequencies  were  8  to  16  kHz. 

2-  The  system  delay  on  the  RM-201  was  set  at  15  millise¬ 
conds,  i.e. ,  the  time  lapse  between  reception  and 
transmission  of  a  signal. 

3.  Velocity  on  the  RM-201  was  set  to  1480  m/s- 

4.  Interrogation  period  was  set  at  10  s. 

5.  Interrogation  frequency  (kHz)  was  set  15  kHz. 

The  numerical  constants  of  (3)  and  (4)  were  used  for  tne 
survey.  Details  of  the  transponders  are  given  ny  Oceano 
Instruments  (1984). 
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II.  THE  OPEB&IIOMAL  PBOCEDOEE  FOR  IMPLANT&TIOH  OF  A 

TRANSPONDER  ARRAY 

A.  SITE  SELECTION 

The  main  criterion  for  the  site  of  the  experiment  was 
that  it  be  relatively  flat  over  an  area  approximately  2  km  x 
2  km.  Since  there  were  no  existing  charts  detailed  enonghr 
a  reconnaissance  survey  was  perfomed  by  the  R/V  ACAMIA  on  10 
and  11  April  1985.  Some  10  north-south  sounding  lines  and 
1-cross  line  were  run  with  1-nmi  spacings  and  -S-min  position 
fixes. 

The  area  chosen  for  the  experiment  (Figure  2-1)  is 
approximately  24  nmi  west  of  Pt-  Lobos,  California,  in  water 
approximately  1600  to  2000  m  deep-  It  is  the  area  from 
36029.5'  to  36031 'N  and  122016*  to  122017. 8'H.  Mini -Ranger 
stations  (Figure  2.2)  were  set  up  at  Point  Sur 
(36018*20. 279"  N  X  121®  53*  56.  179«W)  and  Carmel 
(36033*  49.  176'*  N  X  1  21053*  48. 358'*W)  for  additional  ship 
position  control. 

B.  IllPLANTATION  PROCEDORE 

Tne  transponders  were  deployed  at  pre-selected  sites 
using  a  buoy-first,  anchor-last,  technique.  The  anchors 
used  were  railroad  wheels  which  weigh  about  290  kg  in  water. 
The  mooring  cables  between  the  anchors  and  the  transponders 
and  between  the  transponders  and  the  buoys  are  each  8-m 
long.  The  transponders  weigh  about  46  lbs  in  water,  ana 
each  buoy  has  55  lbs  positive  floatation.  Three  17-in 
Benthos  glass  sphere  floatation  buoys  were  attached  to  each 
transponder  by  means  of  5/32-in  stainless  cables  8-ia  long. 
Altnougn  2  spheres  would  have  provided  sufficent  buoyancy 
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Figure  2-1  The  Project  Area. 
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Figure  2.3  Schematic  of  Bottom  TranspoDder  Mooring. 
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III.  DETEflHIIiAIION  OF  HARMONIC  MEAN  SOOHD  VELOCITY 

A.  INIBODUCTIOH 

Historically,  echo  soundiags  have  been  made  by  assuming 
a  constant,  approximately  average,  value  for  the  velocity  of 
sound  throughout  the  vater  column,  usually 
480  0  ft/s  (=1463  m/s)  or  1500  m/s,  and  then  correcting  for 
variations  of  this  assumed  value  for  the  actual  vater 
column.  Echo  sounders  are  time-measuring  devices.  The 
sounder's  acoustic  transducer  emits  a  sound  pulse  which  is 
reflected  upwards  by  the  bottom- mounted  acoustic  transponder 
and  received  back  at  the  surface  transducer.  The  time  of 
travel  of  the  sound  pulse  is  divided  by  two,  and  this  value 
is  multiplied  by  the  assumed  harmonic  mean  speed  of  sound  in 
sea  water,  thus,  giving  the  distance  according  to  the 
expression  D  =  V  t  .  The  echo  sounder  makes  this  transfor¬ 
mation  electronically  or  mechanically  within  the  device 
itself  and  displays  the  aistance  in  the  water  between  trans¬ 
ducer  and  transponders.  This  distance  is  not  equal  to  the 
true  distance,  since  the  assumed  harmonic  mean  sound  speed 
generally  does  not  equal  the  true  mean  speed  for  a  partic- 
urar  depth.  To  determine  the  true  harmonic  mean  sound  speed 
one  must  know  the  sound  speed  throughout  the  water  column 
and  apply  it  to  solve  for  the  harmonic  mean.  For  a  distance 
measurement  to  be  accurate,  precise  measurements  of  time  and 
the  harmonic  mean  speed  of  sound  in  the  water  column  are 
mandatory. 

B.  THE  VELOCITY  OF  SOOND  IN  THE  OCEAN 

The  velocity  of  sound  in  sea  water  is  a  function  of 
temperature,  pressure,  and  salinity;  V  =  V(T,S,P). 


Because 


tae  sound  velocity  is  not  constant  with  position,  sound  rays 
are  refracted  according  to  Snell’s  law. 

We  have  followed  the  standard  practice  of  the  U.  S. 
national  Ocean  Service  (NOS)  to  use  Wilson’s  {I  960)  equation 
for  tne  speed  of  sound  as  a  function  of  salinity,  tempera¬ 
ture  and  pressure  (Umbach,  1976); 


V-1449 . 14+V^+Vp+Vg+Vg^ 


(3.1) 


wner  e 


V^=4 . 5721T-4 , 4532X10"  V-2 . 6045X10"^T^ 

+7.9851X10"®t'^ 


(3.2) 


Vp»l . 60272X10"^P+1 . 0268X10" V 

+  3 . 5216X10"®P^-3 . 3603X10"^^?"^ 


(3.3) 


Vg=l. 39799(3-35)+!. 69202X10"^(S-35)^ 

V„_-(  3-35 )(-1.1244X10~^T+7. 7711X10"  V 

O  1 « 

+7. 7016X10"®P-1.2943X10^P^+3. 1580X10"®PT 
+  1 . 5790X10"®PT^)+P(-1 . 8607X10"'^T 
+  7 . 4812X10"®T^+4 . 5283X10"®T^ )  ^3  _  5, 

+P^C-2 . 5294X10"^T+1 . 8563X10"°T^ ) 

+P^C-1.9646X10~^°T) . 


la  tiiese  equatioas  the  absolute  pressure,  P,  is  expressed  in 
Jcg/cm2^  temperature,  T,  is  in  °C,  salinity,  S,  is  in  g/kg, 
and  sound  velocity,  V,  is  in  m/s.  For  these  relations  for  V 
the  standard  deviation  from  the  mean  is  0.30  m/s  for  all 
data  obtained  in  the  ranges  -4  <  T  <  30  ®C,  1  <  P  <  1000 
kg/cm2,  and  0  <  S  <  37  g/kg  (Wilson,  1960,  p.  1357). 

C.  HifiHONlC  MEAN  SOOHD  VELOCliy 

The  slant  range  from  a  ship  to  a  transponder  is  deter¬ 
mined  by  measuring  the  transit  time,  t,  of  an  acoustic  pulse 
and  converting  it  to  distance  using  an  appropriate  value  for 
sound  speed.  If  the  mean  sound  velocity,  V,  in  the  water 
column  is  known,  tne  distance,  L,  from  transducer  to  tran¬ 
sponder  can  be  computed  by: 


L  =  (  t/2  )  V 


(3.6) 


m 


For  our  small  project  area  the  distances  between  transpon¬ 
ders  are  less  tuan  2  nmi,  and  it  is  assumed  that  there  is  no 
horizontal  sound  speed  variation,  and  only  corrections  for 
vertical  sound  speed  variation  are  considered.  For  a 
one-way  travel  time  ,  T,  through  a  water  column  of  depth,  Z, 
the  mean  sound  velocity  is; 


V  =  Z  /  T 


(3.7) 


The  meaning  of  "mean  sound  velocity”  is  neither  "mean 
velocity  from  the  surface  to  the  stated  depth”  nor  the  "mean 
value  for  the  velocity  of  sound  through  the  vertical  water 
column".  According  to  Equation  (3.7)  the  mean  velocity  is 
obtained  from  tne  true  depth  divided  by  the  one  way  travel 
time. 
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figure  3.1  A  Layered  Ocean  there  ?;=  (z)  Is  the  Sound 

Speed  in  the  ith  Layer  and  z  Is  Depth. 

consider  nov  a  layered  ocean  having  a  different  sound 
speed  in  each  layer.  Figure  3.1  shows  a  series  of  finite 
layers  of  thickness  az.,  each  with  an  associated  sound 
speed,  Vj^ .  The  time  interval,  AT.,  for  the  sound  wave  to 
pass  vertically  through  the  layer  is  : 

AT.  =  Az.  /  V-  (3.8) 

Summing  tne  time  intervals  for  all  n  layers  in  a  vertical 
water  column  from  surface  transducer  to  bottom  transponder 

T  =  r  at.  =  £  )  (3.9) 

t=l  L=i 

For  a  continuous  function,  v  =  V(z),  this  becomes 


dz/7(z) 


Equation  (3.7)  may  now  be  written  ; 


(3-  10) 


_z _ 

V(2) 


(3.  11) 


Thic  is  the  integral  form  of  the  harmonic  mean-  Using 
Simpson's  Rule,  Equation  (3.10)  can  be  rewritten  as  : 


T  =  (1/V„ +4/V,  +2/72  +4/7,  ♦ 


+  1/V_)  d/3 


(3.  12) 


where  d  =  D/n  (magnitude  of  depth  increments) 

D  =  bottom  depth  at  transponder 
n  =  even  number  of  depth  increments  between 
transducer  and  transponder 
V;  =  measured  sound  speed  at  depth  j. 

Since  the  harmonic  mean  souna  velocity  is  obtained  by 
Equation  (3.11),  tb°  slant  range  from  transducer  to  acoustic 
transponder  can  be  obtained  using  Equations  (3.  10)  and 
(3.11): 

L  =  I  V 

CID  casts  to  140  0  m  were  made  on  18  and  20  May  1  985  in 
tne  project  area  just  prior  to  and  during  our  main  experi¬ 
ment.  Sound  velocities  and  harmonic  means  based  on  the 
sound  velocity  profiles  derived  from  tne  CTD  casts  are 
plotted  in  Figures  3.2  and  3.3  .  Altnough  the  profiles 
Shown  in  these  figures  ace  for  different  days  and  different 
places  within  the  project  area,  they  show  very  little  varia¬ 
tion  in  harmonic  mean  sound  velocity,  and  thus  support  our 
previous  assumption  that  the  sea  water  is  horizontally  homo¬ 
geneous  to  an  acceptable  degree  within  the  study  area- 
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Figure  3,3  Harnonic  Mean  and  Sound  Velocity  Profile  2 


Actually,  the  difference  of  these  tuo  harmonic  mean  veloci¬ 
ties  is  approximately  0.07  m/s  near  the  bottom  as  is  shown 
in  Figure  3.4.  The  mean  of  these  two  harmonic  mean  sound 
velocity  profiles  is  used  in  the  range  calculations  of  the 
following  chapter. 

D.  ACOOSTIC  HATE  fiEFBACTION 

To  simplify  our  acoustic  analysis  it  is  desirable  to 
assume  that  in  the  study  area  sound  follows  straight  paths 
from  the  surface  to  the  ocean  bottom.  To  test  this  approxi¬ 
mation  the  differences  in  path  lengths  from  the  surface  to 
the  bottom  were  calculated  for  varies  depressi-on  angles  at 
the  surface. 

Profiles  such  as  those  of  Figures  3.2  and  3.3  are 
usually  simplified  for  analysis  by  separation  into  an  appro¬ 
priate  number  of  segments  each  having  an  approximately 
constant  gradient  (Figure  3.5).  If  the  velocity  of  sound 
changes  linearly  with  depth,  sound  rays  can  be  shown  to  nave 
a  constant  radius  of  curvature  (Urick,  1983,  p.  124). 
Horizontal  distance,  d,  depth,  z,  radius  of  curvature,  R, 
segments  of  arc,  s,  and  depression  angle.  A,  are  shown  sche¬ 
matically  in  Figure  3.6.  The  path  of  a  ray  when  the  speed 
of  sound  varies  witn  depth  can  be  calculated  by  application 
of  Snell's  law: 

cosAj  /  C,  =  cosA^  /  ^2  “  ***  “  cosAV  (3.13) 

wnere  C<=  sound  speed 

Ai=  depression  angle 

For  an  initial  depression  angle  of  30®  the  profile  of 
Figure  3.5  gives  the  results  shown  in  Table  I. 

The  raaius  of  curvature  can  be  obtained  in  the  following 
manner  (Kinsler,  198  2,  pp399-402).  Let 
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Figure  3.5  Sinplified  Sound  felocity  Profile. 
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TABLE  I 

Depression  Angle  at  Different  Layers 


Points 

1  Sound 

i  Velocity 

1 

1 

Depth 

(m) 

1  Depression 

1  angle  (deg.) 

1. 

1  1493.0 

1 

0. 

1  30. 0 

2. 

i  1483.5 

1 

176.4 

1  30.63 

3. 

1  1484.04 

1 

266.6 

1  30.59 

4. 

1  1481.2 

1 

526.4 

1  30.78 

o. 

1  1432.46 

1 

1097.2 

i  30-69 

6. 

1  1484.34 

1 

1400. 0 

1  30.57 

=  (  c.,,  -  C.  )  /  I  2.^,  -  2.  ) 


(3.14) 


B;  =  -  (C;  /  g.  )  COS  A; 

L  ^  L 


(3.  15) 


wuere 

£•  =  Radius  of  curvature  in  the  depth  increment 

=  Sound  speed  in  the  depth  region 

A;  =  Depression  angle 
g.  =  Speed  gradient 

The  change  in  range  ad.  and  depth  az«  ace  (Kinsler,  1982) 


ad.  =  -E  (  sin  A*  -  sinA^^l  ) 


(3. 16) 


A  2^  =  -E^  (  COSA^^I  -  COSA*  ) 

3  1 


(3.  17) 


so  the  chord  can  be  obtained  by  two  components  of  change  in 
range  and  depth;  that  is. 


D.  =  (  ^d2  ♦  Az.a  )  1/2 

t  Cl 


The  angle,  0;  ,  is 


0;  =  2sin-i  (  D;  /2R;  ) 

*»  L  ^ 


so  the  arc,  s^,  is: 


S  .  =  R;  0; 

<-  L  L 


(3.  18) 


(3.19) 


(3.23) 


where  0*  is  the  angle  in  radians. 

The  range  and  depth  increments  give  tne  straignt  line 


L  =  C  )  2+  (2^z.  )  2]i/: 


The  ray  paths  can  be  obtained  by  summing  arcs; 


S  =  £s; 


(3.21) 


(3.22) 


The  difference  due  to  refraction  between  a  sound  ray's 
actual  path  and  a  straight  line  is  thus: 


E  =  S  -  L 


(3.  23) 


wnere  E  =  the  error  between  two  lines, 

S  =  the  distance  along  the  ray  path,  and 
L  =  the  straight-line  distance  between  the 

ship-mounted  transducer  and  bottom  transponder. 
Ray-trace  computations  were  made  for  various  initial 
depression  angles  (at  the  surface)  down  to  a  depth  of  14 JO  c, 
(Taoles  II  and  III)Since  E<1  cm  for  a  30®  surface 


TABLE  II 

The  Besalt  of  Bay  Trace  Cooputation 


For 

Initial 

Depression  Angle 

=  30  Degrees 

Depth 

Gradient 

Radius  Horiz. 

Depth 

Arc 

Increment 

g 

H  Dist.  (d 

._ImJ _ [mj.  __ 

)  z 

.QnJ _ 

1. 

0.  0538 

32011.36  301-72 

176.4 

349. 5023 

2. 

0.  006  0 

287966.49  152.47 

90.2 

177.  1543 

3. 

0. 0109 

157706-64  437.85 

259.8 

509.  1278 

4. 

0. 0022 

780984-82  960.03 

570.8 

1116. 8993 

5- 

0.0062 

277668-87  511-35 

302.8 

594. 2820 

Total  :  2363.42 

1400.0 

2746.9658 

The 

Distance 

of  Straight  Line 

• 

• 

2746.  9562 

The  Error 

:  0. 

0096  mete 

rs 

depression  angle  and  a  1400-ni  water  depth,  corresponding  to 
a  horizontal  distance  ot  2363  which  covers  the  project 
area,  and  since  the  aepression  angles  encountered  during  the 
project  were  30®  or  greater,  1  cm  represents  an  upper  liait 
to  the  error  resulting  from  our  assumption  of  linear  ray 
paths  from  the  surface  to  the  bottom,  and  a  correction  for 
refraction  is  unnecessary. 


17.  RELATIVE  POSITION  DETER MINATION  OF  ^HE  TRANSPONDER  ARBII 

A.  INTRODUCTION 

Two  methods  to  determine  relative  positions  are 
discussed  in  this  chapter.  They  are  basic  baseline  crossing 
and  the  modified  baseline  crossing  method;  the  former  is 
described  by  Hart  (1  967)  and  Haehnle  (1967),  the  latter  is 
described  by  Fubara  and  Mourad  (1972)  .  The  main  difference 
between  the  methods  is  that  Hart  and  Haehnle  require  accu¬ 
rate  transponder  depths,  but  Fubara  and  Mourad-  do  not. 

The  cloverleaf  metnod  used  to  determine  transponder 
depths  IS  described  by  Haehnle  (1967)  is  introduced  in  tne 
next  section. 

B.  CLOVERLEAF  METHOD  FOR  TRANSPONDER  DEPTH  DETERMINATION 

An  initial  approximation  to  the  depths  of  bottom-mounted 
acoustic  transponders  can  be  obtained  from  a  knowledge  of 
the  water  depth  at  the  point  where  the  mooring  is  released 
from  the  ship.  However,  the  unit  may  move  laterally  while 
descending,  or  bottom  topography  may  have  considerable 
relief  not  evident  on  the  echo  sounder  record.  Hence,  a 
more  accurate  metnod  of  depth  determination  is  required. 

Here  the  depth  is  defined  as  the  distance  between  the  huil- 
mounted  transducer  and  the  bottom-anchored  transponder  wnen 
tne  vessel  is  positioned  exactly  over  the  transponder. 

The  ship  need  not  be  directly  above  the  transponder  but 
should  be  in  close  proximity  to  it.  Although  it  is  very 
difficult  to  cross  over  the  transponder  exactly  with  tne 
ship,  the  depth  may  be  approximated  as  the  minimum  slant 
range  when  the  transponder  is  crossed  over  many  times  using 
a  '‘cloveriedf"  method  (figure  4,1).  Thus,  if  a  transponder 
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The  first  pass  through  the  transponder  location. 

k 


\ 


b.  Maneuver  for  obtaining  four  CPA's. 


Figure  4.1 


Deteraiaation  of  Transponder  Depth. 
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Figure  4.2  Ezanple  of  Depth  Difference. 


is  dt  a  deptu  of  1530  a  and  the  ship  is  displaced  50  meters 
horizontally  from  it,  tne  depth  determined  by  the  cloverleaf 
metncd  would  be  in  error  by  only  0-33  meters  (Figure  4.2). 

The  procedure  used  for  running  a  cloverleaf  maneuver  is 
as  follows  (Haehnle,  1967); 

1-  Each  time  a  transponder  is  dropped  from  the  ship  the 
ship's  position  is  noted. 

2-  The  ship  proceeds  toward  the  transponder  drop  point 
and  passes  near  it;  &  closest  point  of  approacn 
(CPA)  is  then  determined  and  recoraed- 

3-  The  ship  then  turns,  crossing  the  first  tracic  perpen¬ 
dicularly  and  passing  near  the  first  CPA-  Once 
again,  the  CPA  of  tnis  ran  is  cecorded- 

4.  Tne  above  maneuver  is  repeated  until  at  least  four 
such  CPA's  have  oeen  obtained,  each  time  attempting 
to  pass  through  the  previous  CPA  perpendicular  to  the 
tracic. 


Figure  4.J  Saaple  of  Depth ' Determination. 


5.  The  depth  of  the  transponder  is  deterarned  by  plot¬ 
ting  the  slant  ranges  versus  time  for  each  of  tne 
four  CPA  runs.  The  minimum  point  of  tnis  curve  is 
considered  as  the  minimum  slant  range  for  that 
particular  transponder  (Figure  4.3). 

The  S/V  ACA:1IA  ran  similar  cloverleaf  patterns  on  21  Ma 
1985.  Sne  crossed  near  each  transponder  many  times  as 
closely  as  possinie.  The  minimum  distances  and  corrected 
depths  are  listed  in  Table  IV  . 


TABLE  IV 

Transponder  Depth  by  Clorerleaf  Maneuver 
(From  Oceano  £  1985,  p-12]) 


Corrected 

Transponder 

Depth 

(m) 

1 

1366.9 

2 

1331.3 

3 

1381-5 

4 

1406.7 

I _ j 


C.  BASELIBE  CEOSSIBG  MEIHOOS 


Baseline  length  can  he  determined  utilizing  the  baseline 
crossing  method  described  by  Haehnle  (1  967)  and  Hart  (1967), 
which  involves  repeatedly  transiting  the  baseline  while 
simultaneously  measuring  ranges  to  two  transponders.  The 
minimum  sum  of  the  two  horizontal  ranges  as  determined  from 
measured  slant  ranges  and  known  transponder  depths  is  taKen 
to  be  the  baseline  length. 

The  basic  baseline  crossing  method  requires  a  knowledge 
of; 

1.  The  ship's  course  and  speed, 

2.  The  accoustic  slant  range  (or  the  two-way  travel 
times  from  the  ship  to  two  adjoining  ocean- 
bottom  transponders,  and 

3.  The  harmonic  mean  sound  velocity. 


The  ship,  travelling  at  a  constant  speed  and  heading. 


crosses  the  lines  joining  the  transponders  as  shown  in 
Figure  4.4.  The  acoustic  slant  ranges  are  recorded  as  func 
tions  of  time  during  the  crossings.  The  sum  of  two  hori¬ 
zontal  ranges  (the  distance  T1-T2  in  Figure  4.4)  is  minimum 
when  the  ship  is  in  the  vertical  plane  containing  the  two 
transponders.  Similarly,  for  a  f our-transponde r  array,  the 
distances  I1-T2,  12-13,  T3-T4,  T1-T4,  T1-T3  and  T2-T4  are 
determined.  The  minimum  sums  can  be  obtained  by  least- 
sguares  fitting  the  horizontal  ranges  against  time  and 
obtaining  the  vertex  of  the  parabola  as  shown  in  Figure  4.5 
During  the  experiment  four  to  six  crossings  were  made  for 
each  baseline  between  transponders.  The  results  are 
discussed  below  in  Sections  1  and  2. 

Although  the  baseline  crossing  method  is  relatively 
simple  with  respect  to  data  reduction,  it  is  operationally 
time  consuming.  During  the  experiment  it  took  one  hour  and 
thirty  minutes  for  each  baseline  crossing.  Thus  the  eight 
crossings  (4  for  north-south  crossing,  4  for  east-west 
crossing)  took  a  total  of  12  hours. 

1 .  Msic  Baseline  Cr osgin^  Model 

Figure  4.4  shows  a  ship  track  crossing  a  baseline. 
The  ship's  transducer  was  set  to  interrogate  the  bottom 
acoustic  transponders  every  10  seconds  for  this  experiment. 
Sanges  were  obtained  from  each  transponder.  It  is  easily 
seen  (figure  4.4)  that,  as  the  ship  approaches  the  oaseline 
I1T2,  the  total  horzontal  range  at  the  first  ship  point  (0) 
is  greater  than  at  the  second  ship  point  (1)  ,  and  so  on, 
until  a  minimum  sum  is  reached  at  the  instant  of  baseline 
crossing;  from  there  on  the  total  increases  for  each  point 
on  the  track.  Since  the  ship  maintains  a  constant  course 
and  speed,  the  mathematical  function  which  describes  this 
relationship  is  a  parabola  (Figure  4.5)  of  the  form  (Hart, 
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Figure  4.5  Saiple  of  Basge^Time  Parabola  Showing 
Baseline  Crossing  Tine. 


a=AI2  +  BT  +  C 


(4.1) 


wner  e 


A,  B,  C  are  constants  to  be  determined, 

T  AS  ti  me  [  s  ],  and 

P.  is  total  horizontal  range  between  transponders 
£  m]. 

To  determine  the  unknown  variables  A,  3  and  C  in 
Eguation  (4.1)  a  least  squares  adjustment  is  performed.  At 
least  four  sets  of  data  R  =  R(T)  are  needed  to  solve 
Eguation  (4.1).  Since  observed  parameters  contain  random 
error.  Equation  (4.1)  can  be  rewritten  as: 


"N  -V  -V-'f  ."f  w  .lV  .  -l  «  v''-  1%  w  wW  v -j  -J.  -i -i  r-> 


A  T2  ♦  B  T  ♦  C 


H  ♦  V 


(4.2) 


where 


V  is  the  residual  of  the  observation. 
In  matrix  notation  this  corresponds  to: 


fv,' 

f-TA-T,.')' 

-tA'T.,-/ 

/  \ 

A 

Rz 

^  J 

CO  o 

Rri' 

which  may  be  written  as: 


V  + 
UX/ 


T 


R 

nxi 


(4.3) 


(4.4) 


where 


R- 


are  the  residuals  of  the  onservations, 
are  the  numerical  coefficients  of  the 
unknown  variable, 

are  the  unknown  variables (A,  B,  C) ,  and 
are  constant  terms  of  total  horizontal  range. 


The  standard  deviation  of  the  mean  of  sound  velocity 
is  ±0.3  m/s  (Wilson,  19b0).  So,  the  standard  deviation  of 
the  distance  is  the  following,  where  time  is  assumed  error¬ 
less  : 


=  CT^  T  (4.5) 

Each  observation  is  considered  to  be  independent  and 
uncorrelated.  Since  a  measurement  of  high  precision  has  a 
small  variance,  and  one  of  low  precision  has  a  large  vari¬ 
ance,  the  higher  the  weight  the  higher  is  the  precision  and 
vice  versa.  Accordingly,  the  weight  W  of  a  single  obstra- 


tion  Is  defined  as  a  gaantity  that  is  inversely  proportional 
to  the  variances  of  the  observations  {6^)  (Mikhailr  1981, 

p.66)  . 

W  =  d2  /  (5®  (4.6) 

where  6^  =  variance  of  unit  weight 

For  the  adjustment  of  Equation  (4.4),  the  least  squares 
criterion  requires  minimization  of  the  weighted  function  G 
(Mikhail,  1981,  pp.69-73): 

3  =  W,  V2  *  ...  -  (4.7) 

which  in  matrix  form  is  written  as 
t 

G  =  V  W  V  (4.8) 

where 

W 

I  •  >  •  •  Wn.  j 

Now  Equation  (4.4)  can  be  rearranged  as: 

V  =  E  -  T  X  (4.9) 

Substituting  for  V  into  Equation  (4.8),  get 

G  =  (R-TX)^  W(R-TX)  =  (R*-X^T^)  W  (H-TX)  (4.  1J) 

so, 

3  =  R  WR-X  T  WR-E  1»TX  +  X  T^HTX 


^,0  .  .  Q 
0  ^2, 


(4.  11) 


Siace  6  is  a  scalar,  the  eight  hand  side  of  Equation  (4.10) 
is  also  a  scalar.  Furthermore,  the  transpose  of  a  scalar  is 
equal  to  the  scalar: 


k'^^T^WR  =  (x’^T^WR)  =  r'^^MTX 
and 


(4.  12) 


■t 

«  =  W  (4.  13) 

Thus,  Equation  (4.10)  can  be  written  as: 

G  =  F*WR  -  2R^WTX  +  X^(t'^WT)X  ’  (4.14) 


In  Equation  (4.14)  all  matrices  are  constants,  except  X,  the 
matrix  of  unknowns.  For  G  to  ne  a  minimum,  its  partial 
derivative  with  respect  to  X  must  be  equal  to  zero: 


3X 


=  -2R'-  WT  ♦ 


f  t 

2X  (I  WT) 


0 


(4.  15) 


and , 

3xn  3X>x  /txn.  n.x/ 

giving  the  solution  for  X: 

X  =  W  T  )-i  H  R 


(4.  16) 


(4.  17) 


Now,  the  minimum  distance  occurs  when  dR/dT  =  0 


Thus,  to  get  the  aiainani  distance  of  the  sum,  one  may 
substitute  the  time  from  Eguation  (4.  18)  into  Equation 
(4.1)  : 


R  (min.)  =  B2/4A  -  B2/2A  ♦  C  (4.19) 

where  A,  B,  C  are  given  by  Equation  (4.17). 

The  variances  of  A,  B  and  C  (Mikhail,  1981,  p258) 

are : 

6^  =  (T*HT)-»  (4-20) 

The  standara  deviation  of  T  in  Equation  (4.18)  is: 

=  [  (C^^  B/2A2)2  ♦  ((4/2A)2]»/2  (4.21) 

The  standard  deviation  of  the  minimum  distance  E  in  Equation 
(4. 19)  is: 


dL  =[  (T2Cf  )2+(Ti3'  )2+(52+(2AI*3)2d2]i/2  (4.22) 

K  A  o  C  T 

2-  pat  a  Processing 

During  the  experiment  at  sea,  both  north-soutn  anJ 
east-west  baseline  crossings  were  run  (Figures  4.6  and  4.7) 
using  data  collected  by  Oceano  Instruments  (1935a) . 

For  each  baseline  crossing,  11  consecutive  interrogation 
points  were  fit  to  a  parabola  by  least  squares.  The  reason 
for  using  11  points  is  that  this  number  is  sufficient  to 
define  the  parabola  near  its  minimum.  The  minimum  value  of  the 
parabola  can  be  considered  to  define  the  baseline  length- 
The  parabolas  shown  in  the  Figures  4.8,  4.9,  4.10,  4.11, 


Figure  4.6  Morth-Soatii  baseline  Crossing. 

4.12  and  4.13  (at  the  end  of  this  section)  are  for  single 
crossings  between  11-12,  T2-T3,  T3-T4,  T1-T4,  T 1-T3  and 
T2-T4.  Ine  results  are  tabulated  in  detail  in  Table  V,  which 
gives  tne  baseline  lengths  and  their  standard  deviations  for 
each  crossing.  To  get  the  mean  of  these  lengths,  one  can 
use  the  following  eguation; 


h-  =  (El+E2^  •••  ♦Bn)  /  n  =  f  (B) 


(4.23) 


Since  the  standard  deviation  of  the  mean  length  B^in 
Equation  (4.23)  is: 


:L  =[  Idl  df/9hl)2*  {d2  df/dR2)^*  •  •  •  ]  I  /  2 


(4.24) 


we  nave 
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Pigtire  4.7  East-B«st  Baseline  Crossing. 

°m.  =r  <<^Vn)2*(d2/n)  (4.25) 

Equations  4.23  and  4.25  are  used  with  the  data  of 
Table  V  to  compute  the  mean  of  baseline  lengths  and  their 
standard  deviations  (Table  VI)« 

Lyman  BurXe  of  Oceano  Instruments  and  L.  Spielvogel 
of  Seaco,  Inc.,  have  also  computed  the  baseline  lengths. 

They  use  a  different  method  and  technique  using  most  of  the 
data  collected  during  the  experiment.  Their  results  and  the 
differences  between  them  are  tabulated  in  Table  VII-  Due  to 
our  more  limited  data  set  and  uncertainty  in  the  depth  for 
the  baseline  crossing  method,  these  results  differ  from 
those  obtained  with  the  crossing  method  by  about  2  m. 


TABLE  ? 

Baseline  Lengths  Deternined  by  Baseline  Crossings  {a) 


Crossing  J  T1  -  T2  |  T2  -  13  1  T3  -  T4 


1.  I  1793.32±2.43  |  1993.09±0.o1  |  1972-39±2.50 

2-  I  1792-<j2±2.  17  1  1991.92±2.73  |  1969.50+2.02 

3,  (  1793.o7±2.27  {  1992.75±2.38  |  1970.29±2-19 

4.  I  1791.7  7±2.22  |  1994.  33±2.09  |  1969.97±2.23 


Crossing  l  11-14  |  11-13  1  12-14 


1.  I  1860,8  1±3.  25  I  2789.62±2.20  |  2588-66  +  2.03 

2.  I  185a.43±1.81  i  2789.77±1-62  |  2589.70+1.94 

3.  i  1358.63t2.23  1  2789.45±2.  11  |  2587.88±1.14 

4.  I  1858.4  8±2.  08  (  2789.73±1.63  |  2588.  17±2.21 

b.  I  1858.65±2.04  |  2790.22±2.26  |  2589.22±2.2b 

b.  I  !  I  2588-82±2.20 


3-  aodo-f red  Baseline  Crossing  Method 


Mourad  and  Fuhara  (1972)  developed  this  method. 

Their  method  determines  the  relative  positions  of  ocean- 
bottom  transponders  and  the  depth  of  each  transponder. 

There  is  no  necessity  to  accurately  measure  each  depth 
directly;  this  is  the  main  difference  from  the  previous 
work.  The  physical  principle  of  this  method  is  shown  in 
Figure  4.14  S  is  the  baseline  length  and  Li  and  L2  are  a 
pair  of  coplanar  ranges  from  the  ship  to  transponders  1  and 
2  whose  depths  are  Zl  and  Z2,  respectively.  The  mathmatical 
expression  for  this  configuration  is: 


TABLE  ?I 

The  Bean  of  The  Baseline  Lengths 


Mean  Value 
(  m  ) 


Standard  Deviation 
(  m  ) 


T1  - 

T2 

T2  - 

T3 

T3  - 

T4 

T1  - 

T4 

1792. 84 
1992. 42 
1970. 54 
1859- 00 
2789. 76 
2588. 74 


S  =  ( 

L12  - 

Z12  )i 

/  2  ♦ 

(  L22 

Z22  ) 1 / 2 

(4. 

26) 

To  determine 

the 

three 

unJcnowns  S,  Z1  and 

Z2 

by 

least  squares 

method 

,  at 

least 

4 

pairs  of  Li  and 

L2 

are 

needed. 

Now , 

Equati 

on  (4 

.26) 

can 

he  rewritten  as 

f  o 

Hows: 

F  =  S- 

( 

-  Z12 

) 

-  (  h 

22 

-  Z22  )  »/  2 

(4. 

2  7) 

Assume 

A1  =  ( 

L12  - 

Z12  ) 

1/2 

(4. 

28) 

A2  =  ( 

L22  - 

Z22  ) 

1/  2 

(4. 

29) 

Equation 

(4.27)  is  1 

inearized 

and 

partially  difteren 

tinted 
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Figare  4.13  The  Best  Fit  of  Parabola  for  T2-T4  Crossing. 


Figure  4.  14  Baseline  Crossing  Configuration. 

Tae  resultant  observation  equation  for  a  least- squares  solu 
tion  is  (Sourad  ,  1972,  p.21): 

A  A  ♦  B  V  ♦  W  =  0  (4-35) 

which  can  be  written  as: 


Tnus,  the  least- squa  res  solution  (.lourad,  1972, 


p.21)  gives 


A  =-[  P„+A-  (BP-»B'-)-lH 


(4.37) 


V  =-p-iB*  (BP-»B^)-MAA+W) 


(4.33) 


The  variance  of  unit  weight  is^ 


^0  =  ft  (BP-‘B^)-i  (AA+«)  3*M/df}  1/2 


(4.  39) 


The  variance-covariaace  matrix  for  the  adjusted  parameters  X 


£P^  +A*  (BP-ib'‘^)-ia]-i 


(4. 40) 


The  adjusted  parameters  X^.  is; 


Xa  =  X^  +  A 


(4.41) 


where 

A  is  a  coefficient  matrix  for  the  unknowns  correction, 

B  is  a  coefficient  matrix  for  the  observation, 

V  is  a  vector  of  residuals,  representing  the 
corrections  to  observed  ranges,  LI  and  L2 , 

A  is  the  correction  to  assumed  X^  , 

W  is  d  constant  matrix  in  Equation  (4.27),  when 
using  assumed  value  X^  (Sq  ,L1q,L2(,), 

P,  is  the  weignt  matrix  of  L^, 

is  the  weighting  function  associated  witn  X^ , 
dg  is  tne  variance  of  unit  weight 
XoL  is  true  value  of  S,  Z1  and  Z2 ,  and 
X^  is  an  approximate  value  of  Xq,. 

The  weight  matrices  are  taken  to  be  unity  as  the 
measurements  were  of  equal  precision,  lanle  VIII  shows  the 


resulting  horizontal  distances  between  transponders  and  the 
depths  of  each  transponder.  These  results  are  obtained  by 
measuring  pairs  of  slant  ranges  at  the  instant  the  ship 
crosses  the  baseline.  Since  the  ship  crossed  the  baseline 
four  times  for  the  baseline  crossing  method,  there  are  only 
four  pairs  of  slant  ranges  which  can  be  used;  these  repre¬ 
sent  the  smallest  data  sets  needed  for  a  least  squares  solu¬ 
tion. 

The  baseline  lengths  computed  by  Oceano  Instruments 
and  Dr.  Spielvogel  are  given  in  Table  IX  The  difference 
between  Oceano's  and  Spielvogel’s  computed  values  are  witnin 
the  precision  of  the  system.  Wy  results  in  Table  7III 
differ  from  Table IX, since  Oceano  and  Spielvogel  used  much 
larger  data  sets  (more  than  250  data  points) ,  whereas  I  used 
the  much  smaller  sets  (11  points  for  each  of  four  crossings) 
required  by  the  two  baseline  crossing  methods,  namely  eleven 
data  points  for  each  of  four  crossings.  Since  Oceano  and 
Spielvogel  used  a  resection  method  with  four  transponders, 
it  is  not  necessary  for  the  ship  to  travel  the  baseline 
repeatedly:  The  only  requirement  is  that  the  ship  oe  in  the 
vicinity  of  the  center  of  the  experimental  area,  where  data 
can  be  collected  continuously.  On  the  contrary,  the 
crossing  methods  require  a  large  amount  of  ship  time  to 
collect  a  small  amount  of  data.  Thus,  it  is  not  surprising 
that  there  should  be  good  agreement  between  Oceano  and 
Spielvogel  and  rather  poorer  agreement  between  my  results 
and  theirs. 


TABLE  fill 

Hesalts  of  Hodified  Baseline  Crossing  Hethod 


Trans¬ 

ponder 

Pairs 

Horizontal 

Distance  Between 

Transponders  (m> 

Transponder 

Number 

dean 

Depths  (m) 

1-2 

1782-9U9.  51 

1 

1361.  52±4- 53 

2-3 

1993-02±6.30 

2 

1333. 94±2.54 

3-4 

1964.34±8.32 

3 

1375. 76±5.21 

1-4 

1- 3 

2- 4 

1836. 65±9. 21 

2815.47±15.60 

2574.13±7.55 

4 

1407’.  08±3.47 

TABLE  IX 

Conpacison  of  Besalts  for  Different  Hethods 


Lines 

Oceano 

Spielvogel 

Basic 

Baseline 

Crossing 

Modified 

Baseline 

Crossing 

T1-T2 

1792,68 

1791.7 

1792. 84±1.  14 

1782.91±9.51 

T2-T3 

1997.46 

1995.9 

1992.42±0-93 

199  3.02±6.30 

I3-T4 

1973.92 

1972.6 

197  0.54±1.  12 

1964.34±8.32 

T4-T  1 

1661.86 

1860.5 

1859.00±1. 04 

1836.65±9.21 

T1-T3 

2792.69 

2790.8 

2789. 76±0.  89 

281 5.47±15.6 

T2-T4 

2592,18 

2590.3 

2588.74±0.82 

2574. 13±7.55 

V.  CONCLUSIONS  AND  RECOMMENDATIONS 


Two  methods  were  used  to  determine  the  horizontal 
distances  between  the  elements  of  an  array  of  four  acoustic 
transponders  lying  on  the  ocean  bottom  at  an  average  depth 
of  1370  meters,  namely  basic  and  modified  baseline  crossing 
methods.  The  principal  difference  between  the  two  tech¬ 
niques  is  that  the  former  requires  precise  depth  informa¬ 
tion,  whereas  the  latter  does  not.  A  "cloverleaf"  maneuver 
was  described  which  may  be  used  to  determine  the  depths  of 
the  transponders.  Because  a  larger  number  of  acoustic  range 
observations  was  used  for  the  basic  method,  the  standard 
deviations  for  the  horizontal  distances  between  transponders 
using  that  technioue  were  about  one-tenth  those  found  using 
the  modified  method.  Our  results  for  the  basic  method  are 
consistent  with  independent  calculations  made  by  Spielvogel, 
who  used  a  much  larger  data  set.  The  differences  between  • 
our  results  using  the  modified  method  and  Spielvogel 's 
are  within  our  relatively  large  standard  devi ati ons,  wi th  the 
exception  of  the  distance  measured  between  transponders  1 
and  3.  This  discrepancy  is  probably  due  to  the  very  oblique 
crossings  of  the  1-3  baseline  which  were  used. 

In  addition  to  the  smallness  of  the  data  sets  used  for 
our  calculations  other  sources  of  error  in  our  results  may 
stem  from  the  variability  of  the  ship's  speed  and  its  pitch 
and  roll,  which  we  did  not  take  into  account. 

If  the  absolute  positions  of  the  bottom  transponders  are 
required,  the  position  of  the  ship  may  be  determined  by  an 
absolute  positioning  system  such  as  GPS  while  measuring 
slant  ranges  from  the  ship  to  each  transponder.  When  at 
least  three  or  more  absolute  positions  of  the  ship's 
acoustic  transducer  are  known,  the  transponder  coordinates 


may  then  be  adjusted  by  the  method  of  least  squares  to 
obtain  a  best  fit  between  measured  and  calculated  slant 
ranges  from  the  known  transducer  positions.  This  method  is 
introduced  in  the  Appendix. 


4££Ili£II  i 

&BSOLOTE  POSITION  OF  OCEAN  BOTTOM  TBABSPONOEBS 

The  absolate  position  of  a  ship  can  be  obtained  from  the 
Global  Positioning  system  (GPS).  Here,  ship  position  means 
the  position  of  the  GPS  antenna.  The  antenna  position  is  a 
function  of  time.  Once  the  absolate  position  of  the  antenna 
is  obtained,  and  simultaneously,  the  ship's  speed,  heading, 
pitch,  roll  and  acoustic  data,  then  the  absolute  position  ox 
antenna  can  be  used  to  calculate  the  absolute  position  of 
the  bottom  transponders.  The  procedures  are: 

1.  lo  determine  the  absolute  position  of  the  GPS 
antenna: 

2.  To  convert  the  absolute  position  of  antenna  to  the 
hull-mounted  transducer  using  the  ship  heading,  pitch 
and  roll  data;  and 

3.  Simultaneously,  to  convert  the  ship  transducer's 
position  to  the  oottom  transponders  by  using  acoustic 
data. 

A  method  is  introduced  here  to  determine  the  absolute 
position  of  sea-floor  transponders,  if  the  absolute  position 
of  surface  transducer  is  given.  Least-squares  adjustments 
must  also  be  supplied  to  get  the  best  fit  of  the  transponaer 
pos  it  ion. 

The  three-dimensional  geometry  used  is  illustrated  in 

Figure  A.I.  There  are  4  transponders  (Tl,  T2,  T3  and  T4)  . 

Their  estimated  coordinates  are  the  absolute 

position  of  ship-mounted  transducer  is  S  (U>  ,  V-  ,  .rf.  )  ,  and  tne 

J  J  J 

four  slant  ranges  are  1;. 

The  adjusted  distances  I;  between  transducer  and  tran- 
sponders  can  be  computed  given  coordinates  of  any  two 
points: 


mm 


Pigare  A.1  The  fielationship  of  Transducer  (  S  } 
To  Transponders,  (  T-).. 


=  [  (0-Xl)  2+ (  V-Yl)  2*  (H-Z1)  2  (A.  1) 

=  [  (U-X2)  2*  (  V-Y2)  2+ (a-Z2)  2  ]i/a  (A.  2) 

=  [  (U-X3)  2+(V-Y3)2<.  (W’-Z3)  2]i/2  (A. 3) 

%  =  [  (0-X4)  2*  (  V-Y4)  2+ («-z4)  2  ]i/2  (A. 4) 

which  can  be  written  simply  as  follows; 

\  =  [  (Oj-X.  )2*  (V. -y.  )2  +  (w. -Z.  )2]i/2  (A. 5) 


wnere  i  are  1,2, 3, 4.  index  corresponding  to 
each  transponder; 

j  are  1,2  <•••  N.  transducer  positions 
at  the  surface; 

U,V,  W  are  coordinate  of  transducer; 

X,Y,  Z  are  coordinate  of  transponder;  and 

are  adjusted  distance  of  slant  range. 

It  was  necessary  to  determine  the  minimum  number  of 
positions  <Ii)  of  the  hull-mounted  transducer  that  were 
needed  to  fix  the  position  of  four  transponders  in  the 
tar ee-d imensional  coordinate  system.  There  are  three 
unknown  variables  to  be  solved  for  each  transducer  position 
and  there  are  four  transducer  positions  to  solve  for.  To 
solve  for  tnese  twelve  unknown  variables  and  apply  least- 
sguares  adjustment,  toe  numoer  of  equations  must  be  greater 
than  cr  equal  to  the  number  of  unknown  variables.  Tnis  is 
given  by: 


where  4  is  four  eguations  for  each  interrogation; 

N  is  number  of  interrogation;  and 
12  is  twelve  unknown  for  4  transpocders  coordinate. 

The  solution  of  Eguation  (A. 6)  is  N>3.  That  is,  the 
absolute  positions  in  three-dimensional  coordinates  of  all 
transponders  can  be  determined  with  as  few  as  three  consecu 
tive  range  measurements  from  four  fixed  transponders. 

Mikhail  states  the  two-dimensional  distance  condition 
and  its  linearization  (Mikhail,  1981,  pp266-268)  ,  which  can 
be  extented  to  three-dimensional  coordinates.  -  Equation 
(A.  5)  can  oe  linearized: 


l=l0+[  (X-U)  /ip  ]AX*[  (y.-V)  /io  ]AY.+j;  /IpJaZ;  (A.  7) 


The  adjusted  distance  is  then; 

1.  =  ♦  V.  (A.  8) 

where  S-  are  observed  valuesof  the  distance;  and 

V 

v^  are  corresponding  residuals 
Eguation  (A. 7)  then  becomes: 

v-£  (X-M)  /Jo  ]AX-[  (Y-V)  /lo  ]^Y.-[  (Z.-«)  /tp  ]AZ  =  l.o-i.  (A  .  9) 

(.tC  LCU  ULI.'' 

Letting 

f.  =  i.o  -  j?.  (A.  10) 


in  matrix  form  Eguation  (A.  9)  becomes: 


Hi 


K,  b,a  b„  O  0  0  • 

0  0  0  til  O 


4  + 


(A. 11) 


* 

•  J  ^ 


which  may  be  concisely  written  as: 


V  +  B  AX  =  f  , 
4-A/XI  4^^X<2  4^X1 


(A.  12) 


where 

V  is  residual  o£  the  observation; 

B  is  numerical  coefficient  of  the  unknown 
variable  correction  ; 

A  X  is  unknown  variable  correction  ;  and 
f  is  numerical  constant  term  . 

The  weight  matrix  can  be  derived  from  the  standard  devi¬ 
ation  of  the  slant  range,  which  discussed  in  the  previous 
cnapter.  The  least  squares  solution  of  Equation  (A.  12)  is 
(Mikhail,  1981,  p259): 


N  =  B  WB 


t  =  B^Wf 


(A.  13) 


(A.  14) 


AX  =  N-l  t 


(A.  15) 


X  =  Xg  ♦AX 


(A. 16) 


V  =  f  -  BAX 


(A.  17) 


4^ 


WWW 


Wi'i 


(A.  Id) 


82  =  V* 


H  V  /  r 


<A.  19) 


r  =  N  -  3 


(A.  20) 


whero 


N  is  coefficient  matrix  of  the  normal  equations, 

M  is  weight  matrix, 

t  is  vector  of  'constants*  in  the  normal  equations, 

^x  is  vector  of  parameter  corrections, 

Xois  vector  of  approximate  value, 

2  is  vector  of  adjusted  value, 

d^is  standard  deviation  of  adjusted  distance, 

02  is  estimate  of  the  reference  variance, 
r  is  number  of  statical  degrees  of  freedom, 

K  is  number  of  given  transducer  position,  and 

3  is  number  of  observations  necessary  to  specify 

uniquely  the  model  that  underlies  the  adjustment 
problem. 

Since  during  linearization  we  neglect  all  second  and 
higher  order  terms,  we  must  ensure  that  X  is  not  signifi¬ 
cantly  in  error  because  of  this  approximation. 

Consequently,  we  must  iterate  the  solution  by  using  X  as  a 
new  approximation  and  compute  anotner  correction.  This 
procedure  is  repeated  until  the  correction  is  insignificant 
The  final  result  is: 


X.  =  X.  ♦  Ax^=  x„ 
fit  rt*  n,  n 


(A.  21) 


Where 


ms; 


IV*: 


W 


IV' 


n  is  total  iteration, 

is  adjusted  value  at  n  iteration, 

Ax^  is  correction  at  n  iteration  can  be  neglected  compare 
with  previous  correction  or  within  tolerance 
correction,  and 

is  true  value  of  each  variable. 

Since  errors  in  transducer's  position  exist,  it  now 
must  be  established  how  this  error  affects  the  accuracy  of 
transponder  positions.  This  is  discussed  below: 

Let  the  standard  deviations  of  transducer  coordinates  be 


,  ±d^  ,  ±6^.  Recall  Equation  (A.  5)  : 

^^  =  [  (U  -X.  )2*  (V  -Y.  )  2  +  {W  -2.  )  2  ]i/2  (A.  22) 

The  unknown  coordinates  X,Y,Z  ,  are: 

X;^=  U  -  (V-Y.  )  2-(«-z.)  2]»/2  (A. 23) 

Y.  =  V  -  [i2-(W-Z;)  2-(U-X;)2]1/2  (A.  24) 

z.  =  H  -  [l2-(0-X.  )  2-(V-Y.  )  (A.  25) 

^  c  ^  t 

Nov  assume 

FI  =  [  J2-(v-Y.  )  2-(w-Z.)  2]i/2  (A.26) 

F2  =  [I2-(W-Z.  )  2-(U-X.  )  2]i/2  (A.27) 

L  C  U 

F3  =  [i2-(U-X.  )  2-(v-Y.  )  2]i/2  (A. 28) 

u  c 


So,  the  variances  of  Equations  (A.23),  (A. 24)  and  (A, 25)  are 

(Mikhail,  1981,  p.181): 

(V-Y. /F1)2df«-(W-2. /F1)2d2  (A.  29) 

71 


dz  =d2-(j?. /F2i  ad**  (H-Z.  /F2)ada*(a-X. /F2)2d2 

jf  U 


d^a  =d^-(i?^/F3)  a(j^2+  (0-X. /F3)  ad**  (V-Y. /f3)  arfa 
and  the  position  accuracy  is: 
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